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Abstract 

A new approach to CMG (control moment gyro) momentum management and 
attitude control of the space station is developed. The control algorithm utilizes 
both the gravity-gradient and gyroscopic torques to seek torque equilibrium attitude 
in the presence of secular and cyclic disturbances. We show that depending upon 
mission requirements, either pitch attitude or pitch-axis CMG momentum can be 
held constant. We also show that yaw attitude and roll-axis CMG momentum can 
be held constant while roll attitude and yaw-axis CMG momentum cannot be held 
constant. As a result, the overall attitude and CMG momentum oscillations caused 
by cyclic aerodynamic disturbances are minimized. A state feedback controller with 
minimal computer storage requirement for gain scheduling is also developed. The 
overall closed-loop system is stable for ±30% inertia matrix variations and has more 
than ±10 dB and 45° stability margins in each loop. 
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Introduction 


The space station shown in Fig. 1 will employ CMGs (control moment gyros) as a 
primary actuating device during normal flight mode operation. Since the CMGs are 
momentum exchange devices, external control torques must be used to desaturate 
the CMGs, that is, bring the momentum back to nominal value. Some methods for 
unloading CMG momentum include the use of magnetic torques, reaction jets, and 
gravity-gradient torque. For the space station, the gravity-gradient torque approach 
is preferred since it requires no consumables or additional hardware. For this reason, 
various schemes using gravity-gradient torque have been developed [1-9]. 

Two approaches to momentum management using gravity-gradient torque are 
possible. The first of these is a “discrete” or “periodic” momentum control approach 
[1-5] which is basically a feedforward open-loop control scheme. This scheme utilizes 
appropriate samples of CMG momentum, with the sampling frequency of the same 
order as the orbital frequency. This approach, as well as the “predictive” approach 
[7], can handle expected momentum changes; however, it requires accurate inertia 
matrix properties and environmental models to generate the proper attitude steering 
command. In the “discrete” or “periodic” approach, the attitude controller, with a 
much higher control bandwidth, is designed independently. 

The second approach to CMG momentum management is a “continuous” ap- 
proach [6-9] which integrates the momentum management and attitude control 
design. In this continuous, closed-loop control of both the CMG momentum and 
station attitude, the design objective is to establish a proper tradeoff between station 

pointing and CMG momentum management, while satisfying the specific mission 
requirements. 

Expanding on this “continuous” approach, we present a new scheme to the space 
station momentum management and attitude control. The proposed controller pro- 
vides a proper disturbance accommodation. As a result, the cyclic peak of the 
station attitude and CMG momentum oscillation caused by aerodynamic torque is 
minimized. In this paper, CMGs are considered as an ideal torquer; however, the 
CMG gimbal dynamics as well as the CMG steering law [10, 11] should be included 
in the further development of an efficient control law for the overall system. Struc- 
tural flexibility of the space station considered in (12, 13] is neglected because of the 
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low bandwidth nature of the integrated momentum/attitude controller. We assume 
that a strapdown inertial reference system provides relatively noise- free estimates of 
all states and that the body-axis components of CMG momentum can be measured. 
Practical multivariable controller synthesis is accomplished employing various tech- 
niques such as the classical control approaches [12 — 1<], linear-quadratic-regulator 
(LQR) synthesis technique [18, 19], asymptotic disturbance rejection [19, 20], de- 
centralized (partial state feedback) control [21], and robust eigensystem assignment 
techniques [22-24]. A brief summary of a pole placement algorithm for partial state 
feedback controller is given in Appendix. 


Mathematical Models for Controller Synthesis 


The space station in circular orbit is expected to maintain LVLH (local vertical 
and local horizontal) orientation during normal mode operation. The nonlinear 
equations of motion in terms of components along the body-fixed control axes can 
be written as: 


Space Station Dynamics: 
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where 


Ci = — sin 62 cos 0 3 

c 2 = cos 0 \ sin 0 2 sin 03 + sin 0 i cos 0 2 

c 3 = — sin 0 \ sin 0 2 sin 0 3 + cos 0 X cos 02 


Attitude Kinematics (2-3-1 body-axis sequence): 
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CMG Momentum: 
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where (1, 2, 3) are the roll, pitch, and yaw control axes whose origin is fixed at the 
mass center, with the roll axis in the flight direction, the pitch axis perpendicular to 
the orbit plane, and the yaw axis toward the Earth; (0,, 0 2 , 0 3 ) are the roll, pitch, 
yaw Euler angles of the control (body) axes with respect to LVLH axes which rotate 
with the orbital angular velocity, n; (u\, w 2 , u; 3 ) are the body-axis components of the 
absolute angular velocity of the station; (/„, J 22 , I 33 ) are the moments of inertia; I %J 
(i ^ j ) are the products of inertia; (hi, h 2 , h 3 ) are the body-axis components of the 
CMG momentum; (uj, u 2 , u 3 ) are the body-axis components of the control torque 
caused by CMG momentum change; (wj, w 2 , w 3 ) axe the body-axis components of 
the external disturbance torque; and n is the orbital rate of 0.0011 rad/sec. 

For small attitude deviations from LVLH orientation, the linearized equations 
of motion can be written as: 


Space Station Dynamics: 
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Attitude Kinematics: 
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(6a) 

(6b) 

(6c) 



Eliminating (o>i ,u> 2 ,u/ 3 ) from Eqs. 4 and 5, the linearized equations of motion can 
also be written as: 
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Note that the products of inertia cause three-axis coupling as well as a bias 
torque in each axis. Early flight configurations during station assembly may have 
significant misalignment of the control axes with the principal axes. In such cases, 
the above linear, 3-axis coupled equations of motion will be very useful for linear 
control design. It is by now evident that it is very straightforward to derive Eq. 7; 
however, this does not appear to have been done in the literature to the authors’ 
knowledge. 

Fortunately, most practical situations of interest with small products of iner- 
tia permit further simplication in such a way that pitch motion is uncoupled from 
roll/yaw motion. Otherwise, Eq. 7 should be used for 3-axis coupled stability anal- 
ysis and control design. For the case where the control axes are nearly aligned with 
the principal axes (I x = I n , J 2 = J 22 , and I 3 = J33), Eqs. 4 become 


I\Cj\ + n(/ 2 — I 3 )u 3 + 3 n 2 (/ 2 — I 3 )0 x = — u x + w x 

/2W2 + 3 n 2 (I x — J 3 )0 2 = — u 2 4- w 2 (8) 

I 3 u 3 - n(/ 2 - I x )u> x = -u 3 + w 3 


Finally, the linear equations for 0 X ,0 2 , and 0 3 become 

I\9i + 4 n 2 (/ 2 - I 3 )$ x - n(I x - I 2 + I 3 )d 3 = -u x w x (9a) 

I 1 O 2 4- 3 n 2 (I x — I 3 )8 i = — u 2 4- w 2 (9b) 

Iz,8 3 + n 2 (J 2 — /i)^3 4- n (/i — ^2 + I 3 W 1 = ~«3 + w 3 (9c) 

These are the important equations used for the study of passive or active gravity- 
gradient stabilization of Earth-pointing satellites (e.g., see [18, 16]). Since pitch 
motion is uncoupled from roll/yaw motion, pitch control is often treated separately 
from coupled roll/yaw control. 


5 


However, certain configurations of the space station (e.g., the assembly flight 
#3 of Table 1) may need a large pitch TEA because of the small gravity-gradient 
torque available in the pitch axis. In such cases, Eqs. 1 and 2 with small roll/yaw 
attitude errors and small products of inertia become [25. 2G): 

I\0i + (1 + 3 cos 2 0 2 )n 2 (7 2 - I 2 )$ x - n (/, - 7 2 + 7 3 )0 3 

+ 3(/ 2 — 7 3 )n 2 (sin 0 2 cos 9 2 )9 3 = — u, + w x 
I 2 0 2 + 3n 2 (/ l — 7 3 ) sin 9 2 cos 9 2 = — u 2 + w 2 (10) 

^3 + (1 + 3 sin 2 9 2 )n 2 (I 2 - 7 X )0 3 + n(7 x - 7 2 + I 3 )0 X 

d* 3(7 2 — 7i)n 2 (sin0 2 cos^ 2 )0j = — u 3 w 3 

It is evident that roll/ yaw motion is now affected by pitch motion. If pitch attitude 
is held constant with respect to the pitch TEA, the roll/yaw equations can be 
considered as time invariant. Thus, these equations will also be useful for roll/yaw 
controller design for the space station with large pitch TEA. 

Inertia matrices of the Phase I space station as well as the assembly flight #3 
are listed in Table 1. In this paper we are concerned primarily with the phase 1 
configuration. The uncontrolled space station with such inertial properties is in 
an unstable equilibrium when 0, = 0 (i = 1,2,3). Also included are expected 
aerodynamic disturbances which are modeled as bias plus cyclic terms in the body- 
fixed control axes: 

w(<) = Bias + A n sin(n< + <£„) + A 2 „sin(2n< + <£ 2 n) (11) 

The cyclic component at orbital rate is due to the diurnal bulge effect, while the 
cyclic torque at twice the orbital rate is caused by the rotating solar panels. The 
magnitudes and phases of aerodynamic torque in each axis are assumed unknown 
for control design. 


Pitch-Axis Controller Synthesis 

Equation 9b, which is linearized and uncoupled from the roll/yaw equations, is 
used as the basis for pitch control analysis and design. A block diagram of the pitch- 
axis momentum/attitude control loop is shown in Fig. 2. Two control schemes are 


6 


presented here: one without the indicated cyclic -disturbance rejection filter, C 2 (s ), 
and one including the filter. 

Pitch Control without Cyclic-Disturbance Rejection 

Pitch-axis controller consists of a single control input, u 2 and four states, 9 2 , 0 2 , 
h 2 , f h 2 - The pitch control logic is then given by 


u 2 — h 2P & 2 + h 2 D d 2 + K 2 ffh 2 + K 2 [ J h 


(12) 


where the pitch-axis CMG momentum and its integral are included to prevent 
momentum build-up. 

Various methods for the selection of the four gains can be employed. For ex- 
ample, one can place the closed-loop eigenvalues at any desired location, using the 
pole placement technique. The practical problem with this approach is that it is 
not always clear where to place the eigenvalues for satisfactory performance and 
robustness. 


Pitch controller synthesis by iterative successive loop closures is possible, but 
quite tedious. The pitch controller can also be synthesized using an “equivalent” 
compensator. This approach provides physical insight into the proper tradeoff be- 
tween the attitude control and the momentum management. Combining Eqs. 6b 
and 12 to eliminate h 2 gives an equivalent compensator 

s^(K 2 p + sK 2 q) 


u 2 = 


- I< 2H s - K 2I 


9i 


( 13 ) 


Hence, the integrated momentum/ attitude controller can be interpreted as a second- 
order compensator with four parameters. The momentum controller consists of 
double zero at the origin and complex poles. The attitude controller consists of one 
zero on the real axis, which is a conventional proportional-derivative controller. An 
unstable compensator is needed to stabilize the unstable system. It can be easily 
shown that h 2 (s)/w 2 (s) of the closed-loop system has a zero at s = 0. Thus, h 2 (t) 
has a zero steady-state value for a constant disturbance. 

LQR synthesis of the pitch control for feeding back all four states to u 2 can be 
done quickly using a computer code. The LQR technique makes use of a quadratic 
performance index to synthesize state-feedback gains that minimize this index. It 
is especially useful for multivariable systems; all loops are closed simultaneously 
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instead of successively as in classical frequency-domain methods. However, the 
problem with this approach is that the proper selection of the weighting matrices 
is not obvious. It is not always possible to predict the effects of given weighting 
matrices on the closed-loop behavior. 

Hence, a practical approach to the pitch-axis controller design (and the roll/yaw 
controller design of the next section) would be to find the control gains and closed- 
loop poles that result from a wide range of weighting matrices, and simulate the 
corresponding closed-loop system. The gain matrix that produces the closed-loop 
responses satisfying the various requirements becomes the final selection. 

The uncontrolled phase 1 space station has unstable pitch mode with poles at s = 
±1.5n and momentum mode with double pole at s = 0 (one from integral feedback). 
After iterative use of an LQR synthesis code available in the PC-MATLAB and 
CTRL-C software, we select the pitch-axis control gains as listed in Table 2. The 
corresponding closed-loop eigenvalues are: 

-1.6n, -l.On, (—1.5 ± 1.5j)n 

The response of the closed-loop system to the pitch disturbance listed in Table 
1’ with initiai con< iitions of 0 2 ( 0) = 1 deg and 0 2 (O) = 0.001 deg/sec, are shown 
in Fig. 3. The transient responses are satisfactory, while the cyclic aerodynamic 
torque causes the periodic responses of both pitch attitude and pitch-axis CMG 
momentum. The h 2 is bounded with zero mean value, while the d 2 is oscillating 
with respect to 7.5° pitch TEA. The CMG momentum peak and torque demand 
are both well below the allowable limits of about 20,000 ft-lb-sec and 150 ft-lb, 
respectively. 

Pitch Control with Cyclic-Disturbance Rejection 

Depending on the circumstances, either pitch attitude or CMG momentum os- 
cillation, caused by the aerodynamic disturbance torque, may be undesirable. In 

such cases, a cyclic-disturbance rejection filter can be employed as illustrated in 
Fig. 2. 

The proposed cyclic-disturbance rejection filter, C 2 (s), is represented as 

C 2 (s) = n ,v 

* ’ [s 2 + n 2 ][s 2 + (2n) 2 ] ( 14 ) 
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The filter poles will appear in the numerator of the closed-loop transfer function 
02(s)/\v 2 (s) or ^ 2 (-s)/w 2 (s), depending on mode selection. This results in a distur- 
bance rejection at frequencies of n and 2 n for either 0 2 or h 2 . The filter numerator. 
N 2 (s), should be properly designed to stabilize the overall control system. 

For example, the disturbance rejection filter for 0 2 can be represented as 

6l 2 + (n ) 2 « 2 = 0 2 ( 15a) 

02 + (2n) 2 /?2 = $2 (15b) 

where initial conditions for or 2 , 0 2 , d 2 , and 0 2 can be arbitrarily selected (usually, 
zero initial conditions). These filter equations are integrated in the on-board com- 
puter. 

We now repeat the pitch control synthesis of the previous section, augmenting 
the pitch-axis equation with the above filter equations. The pitch control logic, with 
additional disturbance rejection filter states, can then be expressed as 

u 2 = K 2 pB 2 4- K 2 q 0 2 + K 2 f{h 2 + K 2I J h 2 + K 2a oi 2 4- K 2l j,6c 2 4- K 2 S@? + K 2 q02 (16) 

The eight gains of Eq. 16 can be determined using either the pole placement or LQR 
technique. We use the LQR synthesis technique, and the selected control gains and 
the corresponding closed-loop eigenvalues are listed in Tables 3 and 4, respectively. 

Since the disturbance rejection filter is used for pitch attitude, it can be shown 
that ^(•s)/ w 2 ('S) has zeros at s = ±nj and s = ±2nj; hence, an asymptotic dis- 
turbance rejection of a cyclic disturbance of frequencies n and 2n is achieved for 
the pitch attitude. Tra ns fer function h 2 (s)/w 2 (s) also has zeros near s — inj and 
±2nj. This results in minimizing the cyclic peak of the CMG momentum. The 
closed-loop response of this controller is shown in Fig. 4. The CMG momentum is 
bounded with zero mean value, while the pitch attitude is now held constant at 7.5° 
pitch TEA after two orbits. These responses should be compared to those in Fig. 3 
where no disturbance rejection filter was used. It can be seen that the cyclic peak 
of the CMG momentum is also reduced. In other words, the use of disturbance 
rejection filtering results in minimizing the overall oscillations of attitude and CMG 
momentum. 

It is also possible to provide cyclic-disturbance rejection for the pitch-axis CMG 
momentum. In this case, we replace 0 2 of Eqs. 15 by h 2 and then find proper gains of 
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Eq. 16. The closed-loop response for this case is shown in Fig. 5; hence, depending 
on the specific mission requirements, either pitch attitude or CMG momentum can 
be held constant employing the pitch controller illustrated in Fig. 12. 


Roll/Yaw Controller Synthesis 

A design procedure similar to that of the pitch-axis design is followed for the 
roll/yaw controller. Two cases are considered: the first without cyclic-disturbance 
rejection filters, the second with filters. 

Roll/Yaw Control without Cyclic-Disturbance Rejection 

The first case of roll/yaw controller design consists of two inputs, u, and u 3 , 
and eight states, including two integral states for the CMG momentum. The eight 
states are: 

Roll-axis: *1 wj h , J h x 

Yaw-axis: 0 3 u> 3 h 3 J h 3 

Note that we use w x and u ; 3 , instead of 0 X and 9 3 , as state variables. The uncon- 
trolled phase 1 space station with inertial property of Table 1 is unstable in roll/yaw 
and has the open-loop poles of (±1.05 ±0.7 j)n, ±nj, 0, 0. The double pole at s = 0 
are due to the integral feedback of h x and h 3 . Similar to the pitch-axis design, 
the roll/ yaw full-state feedback controller can be designed, iteratively, using the 
LQR or P°l e placement techniques. The multi-input characteristics of the roll/yaw 
axes, however, provide for the calculation of various gain matrices which yield the 
same closed-loop eigenvalues. Although the gain matrix can be completely speci- 
fied by assigning not only the closed-loop eigenvalues but also an allowable set of 
closed-loop eigenvectors, there are many open issues in the eigensystem assignment 
approach [22-24], In our study, we simply use the conventional LQR technique to 
find a closed-loop system satisfying various requirements. 

After iterative use of an LQR synthesis code, we select a 2x8 gain matrix with 
the corresponding closed-loop eigenvalues of (-1.05 ± 0.68 j)n, (-1.04 ± 0.72j')n, 

(-1.42 ± 1.38j)n, and (-1.42 ± 1.38j)n. This choice of closed-loop eigenvalues is 
not claimed to be the best. 
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We now discuss some motivation for developing a decentralized, partial state 
feedback controller for the space station. The configurations of the space station 
will be changing as pavloads/modules are attached or removed. This drives the need 
for a control system design that is either insensitive to the configuration changes 
or able to adapt autonomously. At present, conventional design methods, including 
gain scheduling or mode switching, appear adequate for most of the space station 
flight configurations. Some drawbacks of non-adaptive controllers are that the de- 
sign is often time-consuming, and the controller parameters (e.g., the 1G gains of 
the above roll/yaw controller) must be determined for many operating conditions. 
Furthermore, they must be stored in the on-board computer. 

In order to minimize the computer storage requirement for gain scheduling, 
a roll/yaw control logic with reduced number of gains is considered here. One 
possibility might be to use an LQR synthesis code to find a full-state feedback gain 
matrix and then omit some of the less important gains. However, this approach does 
not always guarantee good performance and stability. Hence, we have developed an 
iterative pole placement algorithm which is an extension of the Bass-Gura formula 
for single-input systems. A brief summary of the algorithm for multi-input systems 
is given in Appendix. 

A decentralized, partial state feedback control logic proposed for the space sta- 
tion is simply: 

«i = + K\dw i + Ki h^i 4- K\i J h x (17a) 

u 3 = + K 3 d<jJ 3 + -K 3///13 -f- K 3I J h 3 (1 7b) 

Notice a particular feedback structure where each control input uses states which 
are more closely related to that control input. We select the eight gains by choosing 
the closed-loop eigenvalues to be the optimal, full-state regulator eigenvalues de- 
termined previously. The eight gains obtained using the pole placement algorithm 
described in Appendix are listed in Table 2. A simulation of the closed-loop system 
with roll/yaw disturbances of Table 1 shows that the transient responses are satis- 
factory and that the steady-state responses are cyclic and remain within specified 
limits. The decentralized controller with eight gains compares quite favorably with 
the full-state feedback controller with sixteen gains. 
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Roll/Yaw Control with Cyclic-Disturbance Rejection 

For the minimization of the steady-state oscillation of roll/yaw attitude and 

CMG momentum, cyclic-disturbance rejection filtering for a multivariable system 
is investigated. 

To clarify some multivariable control issues related to transmission zeros [20], let 
us consider the open-loop transfer function matrix from control inputs to roll/yaw 
attitude and CMG momentum. Combining Eqs. 6a, 6c, 9a, and 9c. we can find the 
transfer function matrix as: 


■ 0,(3) ' 

1 

-[7 3 s 2 + (7 2 - 7,)ji 2 ](s 2 + n 2 ), -[(7, - 7 2 + I 3 )ns]{s 2 + n 2 ) ' 

^3 (■«) 

" A , 

[(A ~ h + 7 3 )ns] (s 2 + n 2 ), -[7,s 2 + 4(7 2 - 7 3 )n 2 ](s 2 + n 2 ) 


Ui(s) - w,(s) 
u 3 (s) - w 3 (s) 


(IS) 


hi(s) 

_ 1 

s, n 


u,(s) 

h 3 (s) J 

s 2 + n 2 

. 5 J 


. u a00 . 


(19) 


where A = 7,7 3 (s 2 + n 2 )[s 4 + (1 + 3*:, + M 3 )nV + 4*,* 3 n 4 ], * (f 2 - 7 3 )/7,, fc 3 = 

(h-h )// 3 - 

It is apparent that ±nj are transmission zeros of Eq. 18. We also notice from 
Eq. 19 that resonance of CMG momentum can happen for sinusoidal control inputs 
of frequency n. In other words, a cyclic-disturbance at orbital rate cannot be 
rejected for both roll and yaw attitude. However, in this case, where we use the 
CMG momentum dynamics represented as Eqs. 6, it is not evident why a cyclic- 
disturbance at orbital rate can be rejected for the yaw atttiude and not for the 
roll attitude. In order to prove such inherent characteristics of roll/yaw dynamics, 
the following CMG momentum dynamics with a proper modification of Eqs. 9 are 
investigated: 


hi = u, 

h 3 = ti 3 


For this case, the following transfer function matrix can be obtained: 



1 

. . 

~ A 


-7 3 s 2 (s 2 + n 2 ), 
ns(7 3 - 7 2 )[s 2 + (2n) 2 ], 


(7 2 - 7,)ns(s 2 + n 2 ) 
— s 2 [7is 2 + (7, + 37 2 — 37 3 )n 2 ] 


(20a) 

(20b) 
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( 21 ) 


“i(-s) 
u 3 (s) _ 


' / 3 s 2 + (/ 2 -/,)n 2 . 

(I\ — I 2 4- /i)ns 

w,(.s) 

— (/l “ I 2 + I^TIS, 

I\S 2 + 4(12 — ^3 

. W 3< 5 ) . 


M 5 ) 

1 

‘ 1, 0 ‘ 

tii(s) 

. h 3 (s) 

s 

. °, 1 . 

. u ^( s ) 


where A = IJ 3 s 2 [s 4 + (1 + 3^ + k { k 3 )n 2 s 2 + = (I\ - / 3 )//,.Jt 3 = 

{I 2 -h)fh. 

It can be shown that ±nj are transmission zeros of the transfer function matrix 
from (u x ,u 3 ) to {O\,0 3 ). The zero at s = ±nj appears in both the ^i(s)/«j(s) 
and ^i(^)/u 3 (s) transfer functions, but not in yaw attitude channels. Therefore, a 
cyclic-disturbance rejection at the orbital rate is not possible for roll attitude, while 
it is possible for yaw attitude. Furthermore, it can be shown that ±nj are not 
transmission zeros of the transfer function matrix from (u x ,u 3 ) to (h\,9 3 ). Based 
on these inherent properties of the roll/yaw dynamics, a cyclic-disturbance rejection 
for the roll-axis CMG momentum and yaw attitude is considered. 

The roll/yaw cyclic-disturbance rejection filters proposed for the space station 

are: 


Qi + ( n ) 2 a x = /i x (23a) 

0i + (2n) 2 /?i = h x (23b) 

a 3 + ( n ) 2 or 3 = 6 3 (23c) 

03 + (2n) 2 /? 3 = $ 3 (23d) 


Using a LQR synthesis code, we can determine a 2x16 gain matrix for the roll/yaw 
equations augmented by the above filter equations. However, we consider here a 
decentralized, partial state feedback controller in order to minimize the computer 
storage requirement for gain scheduling. 

A decentralized control logic with cyclic-disturbance accommodation is given 
by: 


u x = I\\pO\ + K\duj\ + Kifjhi + K 11 J h\ 

+ A Xa ar x + AiqQi + K\p0x + K l p0\ (24a ) 

u 3 = K 3 p6 3 + K 3 d^>3 + K 3 [fh 3 + K 3I J h 3 

+ Kza° t 3 + A' 3 ad 3 + K 3 00 3 + K 3 q0 3 (24b) 
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Proper closed-loop eigenvalues are first selected from LQR design iterations, and the 
sixteen gams of Eqs. 24 are then determined by using the pole placement algorithm 
derived in Appendix. These gains and closed-loop eigenvalues are listed in Tables 3 
and 4, respectively. It is emphasized that the decentralized control is not necessary 
if a state feedback controller with 32 gains can be easily implemented without any 
on-board computer storage problem for many different operating conditions. The 
decentralized controller is proposed here as an additional capability of our control 
system for the space station. 

Bode plots of the closed-loop transfer functions are shown in Fig. 6. Both 
*iO)/w i(s) and 0 3 (s)/\v 3 (s) have zeros at s = ±nj and ±2 nj; hence, an asymptotic 
disturbance rejection for the roll-axis CMG momentum and yaw attitude is achieved 
for cyclic disturbances of frequencies n and 2 n. It can be shown that ±nj are 
blocking zeros of the closed-loop transfer function matrix from (w,, w 3 ) to (h,,8 3 ). 
Also 0,(s)/vf x (s) and h 3 (s)/w 3 (s) have zeros near s = ±nj and ±2n>, which results 
in minimizing the cyclic peak of the roll attitude and yaw-axis CMG momentum. 
The responses to the roll/yaw disturbances listed in Table 1, with initial conditions 
of 9,(0) = 0 3 (0) = 1 deg and u/^O) = u 3 (0) = 0.001 deg/sec, are shown in Fig. 7. 
The roll-axis CMG momentum approaches the zero steady-state value, while roll 
attitude oscillates at orbital rate. The yaw attitude approches a constant steady- 
state value, while the yaw-axis CMG momentum oscillates at orbital rate. As a 
result, the overall attitude and CMG momentum oscillations are minimized. In fact, 
the responses of the decentralized controller with 16 gains compare quite favorably 
with those of a full-state feedback controller with 32 gains. 

From simulations, it was verified that the overall closed-loop system, with the 
feedback gains of Table 3, is stable for ±30% inertia matrix variations, and that it 
has more than ±10 dB and 45° stability margins for each loop. 


Conclusions 


We have presented a new approach to the space station momentum manage- 
ment and attitude control. The proposed control algorithm provides a proper 
cyclic-disturbance accommodation. We have shown that either pitch attitude or 
pitch-axis CMG momentum, depending on the specific mission requirements, can 
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be held constant. We have also found some inherent physical property of the cou- 
pled roll/yaw dynamics in terms of a transmission zero of a multivariable svstem. A 
decentralized, partial state feedback controller which requires minimum computer 
storage requirement for gain scheduling was also developed. 


Appendix: Pole Placement Algorithm for Partial State Feedback 

A pole placement algorithm for the decentralized, partial state feedback control 
design is presented here. This iterative algorithm for a multi-input system is an 
extension of the Bass-Gura formula for a single-input system. 

Consider a linear, time-invariant system with a single input described by: 

x = Ax + bu 
u = — kx 


where x is the n-dimensional state vector, u is the scalar input, A is an nxn matrix, 
b is an nxl column vector, and k is a lxn row vector. 

The Bass-Gura formula for determining the feedback gain, k, of a controllable 
system is then given by 


a* - a<> = kCT 


where 


A(s) = characteristic polynomial = s n + ajs " -1 -| + a n _,s + a n 

a<> = [ a io, 02 0 , • ■ • , a no ] for the open-loop system 
a c = [a lc , a 2c , • ■ • , a nc ] for the closed-loop system 
C = controllability matrix = [b, Ab, • • • , A n- 1 b] 

T = upper-Toeplitz matrix whose first row is [ 1 , cq, • • • , a n _[] 

Modification can be made to the Bass-Gura formula for a multi-input system 
described by: 


x = Ax + Bu 
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u = - Kx 


where x is the n-dimensional state vector, u is the in-dimensional input vector. A 
is an nxn matrix, B is an nxm matrix, and K is a mxn gam matrix. 

Let b, be the i-th column of B; that is, B = [bi, bj, • • • , b m ]. Also let k, be the 
i-th row of K. Defining be the closed-loop system matrix at 7-th iteration, the 
algorithm can now be described as 

x = A (j) x + £b,dui J+,) 

i 

+ du\^ 

where du\ J+1 ^ is used in the (7 -f l)-th iteration to improve the convergence to the 
desired eigenvalues. The incremental feedback gain is related to dup* as 

du'p = — dk^x 
k- j) = k- J_1) + dkj - '* 

Small changes in the feedback gains are related to the small changes in the charac- 
teristic polynomial coefficients as 

da (j) = [dk ( 1 j) C ( i-' ) H l-dk^C^] T (j) 

where 

C< J) = [b,-, A< J >b„ ■ • • , ( A^) n - 1 b,] (t = 1, • • • , m) 

T <j) = upper-Toeplitz matrix with the first-row, [ 1 , 0 ^, • • • , 

The above iterative algorithm can be considered as a general form for the pole 
placement for multi-input systems. For the full-state feedback problem, the choice 
of the gains are not unique. There are some methods for the unique selection of 
the gains, such as eigensystem assignment approaches (e.g, see [22-24]). In this 
paper, we consider the case where the number of total feedback gains are specified 
by the number of the states. Then a ps controllability matrix is formed from 
the selected rows of C (j) matrices by replacing the unselected gain components with 
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zeros. This approach is only valid when the system is controllable by the selected 
gain components. 
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Table 1. Space station parameters 


Parameters 

Assembly 

Flight#3 

Phase 1 

Inertia (slug-ft 2 ) 


Ai 

23.22E6 

50.28E6 

I 22 

1.30E6 

10.80E6 

^33 

23.23E6 

58.57E6 

In 

— 0.023E6 

— 0.39E6 

^13 

0.477E6 

A0T6E6 1 

I 23 

— 0.011E6 

0.16E6 

Aerodynamic torque (ft-lb) for Phase 1 

w i 

1 + sin( nt) + 0.5sin(2n<) 

w 2 

4 + 2 sin (nt) + 0.5 sin(2 nt) 

W 3 

1 + sin (nt) + 0.5sin(2ni) 


Table 2. Controller gains without cyclic-disturbance 
rejection filter for the Phase I space station 


Gain 

Pitch 

Roll 

Yaw 

Kp (ft-lb/rad) 

2.448E+2 

1.113E+3 

1.173E+3 

Kd (ft-lb/rad/sec) 

1.465E+5 

5.161E+5 

1.283E+6 

Kh (ft-lb/ft-lb-sec) 

7.523E— 3 

5.935E— 3 

1.540E— 2 

Kj (ft-lb/ft-lb-sec 2 ) 

3.546E— 6 

1.105E— 5 

1.990E— 6 
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Table 3. Controller gains with cyclic-disturbance rejection 
filter for the Phase I space station 


Gain 

Pitch 

Roll 

Yaw 

Kp (ft-lb/rad) 

3.425E+2 

7.026E+2 

9.254E+2 

A’d (ft-lb/rad/sec) 

1.972E+5 

4.382E+5 

4.166E+5 

Kh (ft-lb/ft-lb-sec) 

1.0S9E— 2 

4.779E— 3 

1.817E— 4 

Ki (ft-lb/ft-lb-sec 2 ) 

3.953E— 6 

— 8.170E— 7 

— 1.469E— 6 

K a (ft-lb/rad) 

1.006E— 4 

4.383E— 10* 

— 3.375E— 4 

Ki, (ft-lb/rad/sec) 

5.659E— 2 

7.904E— 7* 

6.606E— 3 •< 

K g (ft-lb/rad) 

1.434E-5 

— 3.309E— 9* 

— 4.503E — 4 

Kp (ft-lb/rad/sec) 

7.608E— 2 

— 4.940E— 7* 

2.412E— 1 


* in units of ft-lb/ft-lb-sec 3 , ft-lb/ft-lb-sec 2 , ft-lb/ft-lb-sec 3 , 
and ft-lb/ft-lb-sec 2 , respectively. 


Table 4. Closed-loop eigenvalues for the Phase I in units of orbital rate, 
n = 0.0011 rad/sec 



Momentum/ Attitude 



Disturbance Filter 

Pitch 


-0.3 ± 1.0 j 

-0.3 ± 2.0 j 

Roll /Yaw 

-0.23, -0.68 -0.66 ± 1.51 j 

-0.23 ± 0.92 j 

—0.20 ± 2.02/ 


— 1.02 ± 0.29/ -1.50 ±0.84 j 

—0.26 ± 1.04j 

-0.62 ± 2.29; 
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Figure Captions 


Figuer 1 The Phase I space Station 

Figure 2 The pitch-axis momentum/attitude control system 

Figure 3 Pitch-axis response without disturbance rejection filterin 

Figure 4 Pitch-axis response with disturbance rejection filtering 
for the pitch attitude 

Figure 5 Pitch-axis response with disturbance rejection filtering 
for the pitch-axis CMG momentum 

Figure 6 Closed-loop frequency response for the roll /yaw axes with 
disturbance rejection filtering for the yaw attitude and 
roll -axis CMG momentum 

Figure 7 Rol 1 and yaw 
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